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Abstract 

The virtual photon structure function gJ(x,Q 2 , P 2 ), which can be 
obtained in polarized e + e~ colliding-beam experiments, is investigated 
for A 2 < P 2 C Q 2 , where — Q 2 (— P 2 ) is the mass squared of the probe 
(target) photon. The analysis is made to next-to-leading order in QCD, 
in the framework of the QCD improved parton model with the DGLAP 
evolution equations. The non-leading corrections significantly modify 
the leading log result, in particular, at large x as well as at small x. 
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In the last several years, the nulcleon spin structure functions g^ n \x, Q 2 ) and 
9% \ x , Q 2 ) have been extensively studied by deep-inelastic scattering of polarized 
leptons on polarized nucleon targets. The information on the spin structure of pho- 
ton would be provided by the resolved-photon process in polarized electron and pro- 
ton collision in the polarized version of HERA|I[. More directly, the spin-dependent 
structure function of photon can be measured by the polarized e + e~ collision in the 
future linear colliders. For real photon (P 2 = 0) target, there exists only one spin 
struture function, gj (x, Q 2 ), which is equivalent to the structure function W^(x, Q 2 ) 
(gl = 2W2) discussed some time ago in @. The leading order (LO) QCD correction 
to g\ was first studied in || and later in @, ||. Recently the first moment of gj 
attracted attention in the literatures ||, |6], [7] in connection with its relevance to 
the axial anomaly, which also plays an important role in the QCD analysis of the 
nulcleon spin structure function . 

Now the next-to-leading order (NLO) QCD analysis for gl is possible, since the 
required spin-dependent two-loop splitting functions in DGLAP evolution equations 
or, equivalently, the two-loop anomalous dimensions of the relevant operators have 
been calculated recently || |9| . Actually there has been already an analysis of gl for 
the case of the real photon target by Stratmann and Vogelsang ||10|| . 

In this paper we shall investigate the polarized virtual photon structure function 
gj(x, Q 2 , P 2 ) to the NLO in QCD, where — Q 2 (-P 2 ) is the mass squared of the 
probe (target) photon (Fig.l). We examine gj (x, Q 2 , P 2 ) for the kinematical region: 
A 2 < P 2 < Q 2 , where A is the QCD scale parameter. The advantage in studying 
the virtual photon target is that we can calculate the whole structure function 
entirely up to NLO by the perturbative method. On the other hand, for the real 
photon target there appear the non-perturbative pieces in NLO [jllj, which we may 
dispose of by using, say, vector-meson-dominance model. The unpolarized virtual 
photon structure functions F% (x, Q 2 , P 2 ), and F2{x, Q 2 , P 2 ) were studied in LO 
and in NLO |13| . 

The analysis can be performed either in the framework of the QCD improved par- 
ton model |14| using the DGLAP parton evolution equations or in the framework of 
the operator product expansion (OPE) supplemented by the renormalization group 
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(RG) method. In this paper we follow the former approach. 

Let g^(a;,Q 2 ,P 2 ), GF±(x, Q 2 , P 2 ), Tl(x,Q 2 ,P 2 ) be quark with z-flavor, gluon, 
and photon distribution functions with ± helicities of the longitudinally polarized 
virtual photon with mass — P 2 . Then the spin- dependent parton distributions are 
defined as 

Aq* = q\ + q\ - qi - qt (1) 
AG 7 = G\ - Gl , Ar 7 = r 7 - n . (2) 

In the leading order of the electromagnetic coupling constant, a = e 2 /47r, Ar 7 
does not evolve with Q 2 and is set to be Ar 7 (x, Q 2 , P 2 ) = 5(1 — x). For later 
convenience we use, instead of Aq\ the flavor singlet and non-singlet combinations 
of spin-dependent quark distributions as follows: 

Ag 7 = £ Ag* (3) 

i 

Ag 7 s = £e?(V-^2), (4 ) 

where rif is the number of flavors of active quarks. In terms of these parton dis- 
tribution functions, the polarized virtual photon structure function gj(x, Q 2 , P 2 ) is 
expressed in the QCD improved parton model as[] 



9l{x, Q\ P 2 ) = f ~ UlUv, Q\ P 2 ) C S (~, Q 2 ) + AG\y, Q 2 , P 2 ) C G (*, Q 2 ) 

J* y l v y 

+Aq NS (y,Q 2 ,P 2 ) C NS (-,Q 2 )\ + C^(x,Q 2 ) , (5) 

y J 

where C S (C NS ), C G , and G 7 are the coefficient functions corresponding to singlet(non- 
singlet)-quark, gluon, and photon, respectively, and are independent of P 2 . 

The parton distributions Aq$, AG 7 , and AqJj S evolve with Q 2 . Their evolution 
equations are expressed in a compact matrix form, once we introduce a row vector 
Aq 7 = (Ag 7 , AG 7 , Aql s ): 

dA * { *'%' P2) = A*(*. Q 2 ) + f 1 ^ Aq 7 (y, Q 2 , P 2 ) AP (^ Q 2 ) , (6) 
a mQ z Jx y y 



1 The detailed derivation of this formula will be reported elsewhere flf 
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where the elements of a row vector Ak = (AK$, AK G , AKns) are polarized split- 
ting functions from virtual-photon to singlet-quark, gluon, and non-singlet-quark, 
respectively. The 3x3 matrix AP(z, Q 2 ) is written as 

(AP q s q (z,Q 2 ) AP Gq (z,Q 2 ) \ 

AP(z,Q 2 )= AP qG (z,Q 2 ) AP GG (z,Q 2 ) . (7) 

V AP q N q s (z,Q 2 )) 

with its element APab representing a splitting function of 5-parton to A-parton. 

Once we get the information on the coefficient functions in Eq.flU) and parton 
splitting functions in Eq.(||), we can predict the behaviour of gj(x, Q 2 , P 2 ) in QCD. 
The NLO analysis is now possible since the spin-dependent one-loop coefficient func- 
tions and two- loop parton splitting functions are available [§, |9|. There are two 
methods to obtain gj(x, Q 2 , P 2 ) in NLO. In one method, we use the parton splitting 
functions up to two-loop level and we solve numerically Aqr 7 (a;, Q 2 , P 2 ) in Eq.(|6]) 
by iteration, starting from the initial quark and gluon distributions in the virtual 
photon at Q 2 = P 2 . The interesting point of studying the virtual photon with mass 
— P 2 is that when P 2 ^> A 2 , the initial parton distributions of the photon are com- 
pletely known up to NLO in QCD. Then inserting the solved Aq 1 (x, Q 2 , P 2 ) into 
Eq.(|5|), and together with the known one- loop coefficient functions we can predict 
gJ{x,Q 2 ,P 2 ) in NLO. 

The other method, which is more common than the former, is by making use of 
the inverse Mellin transformation. We follow the latter method in this paper. First 
we take the Mellin moments of Eq.(|5|), 

f 1 dxx n - l g!{x, Q 2 , P 2 ) = Aqi(n, Q 2 , P 2 ) ■ C(n, Q 2 ) + CI (8) 
Jo 

where we have defined the moments of an arbitrary function f(x) as 

f(n)= C dxx n - x f{x). (9) 



Henceforce we omit the obvious n-dependence for simplicity. We expand the splitting 
functions Ak(Q 2 ) and AP(Q 2 ) in powers of the QCD and QED coupling constants 
as 

Ak(Q 2 ) = >(°» + ^Afc« + - (10) 
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AP(Q 2 



a s {Q 2 
2tt 



« S (Q 2 
2tt 



AP« + ■ 



(11) 



and introduce t instead of Q 2 as the evolution variable [16 



Aln^. 

A « s (g 2 ) 



(12) 



Then, taking the Mellin moments of the both sides in Eq.(|6|), we find that Aq 7 (t)(= 
Aqr 7 (n, Q 2 , P 2 )) satisfies the following inhomogenious differential equation p 



dAq\t) 
Jt 



a 2tt 
2n\ a s 



2/3 



J AP (o) 



Of. 



+Aq^(t)< AP W + n 

[ 27T 



+ 0(a 2 s 



(13) 



where we have used the fact the QCD effective coupling constant a s (Q 2 ) satisfies 



da s (Q 2 
dlnQ 2 



-Pi 



a s (Q 



2\2 



4tt 



(14) 



with /3 = 11 - 2n//3 and ft = 102 - 38n//3. Note that the P 2 dependence of Aq 7 
solely comes from the initial condition (or boundary condition) as we will see below. 
We look for the solution Aqr 7 (t) in the following form: 



Aq\t) = Aq<^\t) + Aq^(t) 



7(1) ( 



(15) 



where the first (second) term represents the solution in the LO (NLO). First we 
discuss about the initial condition of Aq 1 . In the framework of OPE, the twist- 
2 operators R l n (i = S, G, NS, 7) are relevant for the polarized photon structure 
function g[. The expressions of R l n are found, for example, in Ref.||. For — p 2 = 
P 2 ^> A 2 we can calculate perturbatively the photon matrix elements of the hadronic 
operators R l n (i = S,G,NS). Choosing the renormalization point /i 2 to be close to 



P , we get, to the lowest order [13 



( 7 (P) I KM I 1(P)) = ^{-\ K °n + A n)> * = S,G,NS 



(16) 
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where K®' 1 = (K^) 1 are one-loop anomalous dimension matrices between the photon 
and hadronic operators. In fact the K^-terms and A^-terms represent the opera- 
tor mixing between the photon operators and hadronic operators in LO and NLO, 
respectively, and the operator mixing implies that there exists quark and gluon dis- 
tributions in the photon. When we renormalize the photon matrix elements of the 
hadronic operators at fi 2 = P 2 , we obtain 

a 



{i{ P )\K{n)\i{v))\ i 



,2_pa 



A 1 



(17) 



which shows that, at fi 2 = P 2 , quark distribution exists in the photon, not in the 
LO but in the NLO. Thus we have 



Aq 7(0) (0) = 0, 
Explicitly, A n is given by 



a 



Aq^(0) = -A n . 



6 (< e l >, 0, < e 4 > 



< e 2 > 2 ) Air 



V? (19) 

where < e 2 >= Y.i e V n f an d < e4 >= J2i e i/ n f, an d A™ G is the finite part of the 
gluon matrix element of the flavor-singlet quark operator, whose expression in the 
MS scheme is given by ||, 

4 1 1 

n . 



A n 



In 



n 



n(n + 



1 " 1 



n + 1) 



n- 



(20) 



With these initial conditions, we obtain for the solution Aq"'(t) of Eq. (|T3D to the 
NLO (see fT§ for details): 



Aq\t)/ 



a 



^ TS-0 \ ~> pn 1 



1 - 



y p n — - — + — k° y p n (i 



Us(t) 

a s (0) 
1 



i+ 



A?/2/? 



P" P n 1 
— K° J 1n _J 1 

n ^2#, + A?-A? A?/2A) 



a.(0) 



2^ 



+ 2(3 A n . 



pn ^(1) pn I 

2, '71 7 



A \^ nn A"/2/3( 



2/3 + a; 



5 ' 2/3 ^ u 



n / j i 







1 — 

1 ^ Wo 



a s (t) 



"2/3 



(21) 
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Here the moments of the splitting functions are related to the anomalous dimensions 
of operators as follows [§. [|: 



AP (o) 
Afc<°> 



AP (i) 



1 



In 



4 «' 



1 



and we have introduced the projection operators given by 



AP (o) 



pn pn 

i 3 



-i E 



+ : 



i^j, 
P? i = j ' 



E^ r 



(22) 
(23) 

(24) 
(25) 



where A" are the eigenvalues of the matrix 7°. In fact, the parton distributions 
A<jr 7 (t) , do depend on the initial conditions Aq 7 (0) = (a/4Tr)A n but the structure 
function gJ(x,Q 2 , P 2 ) itself is independent of Aq 7 (0) in NLO in QCD, which will 
be discussed in detail in 

When the coefficient functions are given, up to one-loop level, by, 

a s (Q 2 ) 



C s (n,Q 2 ) 

C G (n,Q 2 ) 

C NS (n } Q 2 ) 
C y (n,Q 2 ) 



< e z > 1 



< e 2 > {0 + 



47T 

«5(Q 2 



47T 



-52 



47T 



B 



NS 



a 



47T 3 7 



(26) 

(27) 

(28) 
(29) 



then we obtain from Eq.(^) the following formula for the moments of gj(x, Q 2 , P 2 ) 
in the NLO: 



[ l dxx n ~ l gUx,Q 2 : P 2 ) 
Jo 

V P n 1 47r 

^ 4 l + Ar/2/3 a s (Q 2 



a 1 



4tt 2/3, 



i=+,-,NS 

+ E ^ 

i=+,-,NS 

+ E 

i=+,-,NS 



MP" 



A-/2/3 +l ' 



^(Q 2 ) 
,« S (P 2 ) 

a s (Q 2 ) ' 



A^/2/3 ' 



A"/2/3 +l ' 



+ C n + C(a s 



(30) 
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where P™, Af, B™ and C n are given by 
Pf = K° n I?C n Q.,0) 

pn * (1) pn 



^ 2p + A™ - A* 

C n (l, 0)-^— - 2A>A,/f C B (1, 0) 



1 _ K o^_pn r n m 1 - Ar/2/3o 



(31) 



(32) 



pn^(l) pn ■) 



C n = 2[3 {3n f <e*>B^ + A n -C n (l,0)} 



+ K°P n 



( <e 2 >B n s 
<e 2 > B% 



1 



1 + A"/2/3 



(33) 



(34) 



and the column vector C n (l, 0) = (< e 2 >,0,1) T represents tree-level coefficient 
functions. The expressions of Eqs. (|30|) and ([3~T|)-(j34T) are actually the same in form 
as the ones obtained before by one of the authors and Walsh for the case of the 



virtual photon structure function F% [|13]. All the quantities necessary to evaluate 



P™, A™, B™, and C n are now known, which will be reported in [15]. Especially we 
have at hand the results of the two-loop anomalous dimensions || |9| and one-loop 
coefficient functions Bf (I = S, G, NS, j) in Eqs.(||)-(||) 0, |, |] which were 
calculated in the MS scheme. 

In the case of real photon target, the polarized structure function gJ(x,Q 2 ) 
satisfies a remarkable sum rule which holds true in all order of a and as [7[ , 



( l gJ(x,Q 2 )dx = 0. 
Jo 



(35) 



Now we can ask what happens to the first moment of the virtual photon structure 
function g\ (x , Q 2 , P 2 ) . This can be studied by taking the n — > 1 limit of Eq. fl30|) . 
Note that we have the following eigenvalues for the one-loop anomalous dimension 



matrix 7,' 



n=l- 



A 



n=l 



0, A 



n=l 



-2/3 , 



\ n=l 
A NS 







(36) 
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Physically speaking, the zero eigenvalues A™ =1 = A^. 1 = correspond to the 
conservation of the axial-vector current at one-loop order. The other eigenvalue 
A™ =1 = — 2/3o, which is negative, is rather an artifact of continuation of the gluon 
anomalous dimension to n = 1, since no gauge- invariant twist-2 gluon operator ex- 
ists for n — 1. However, in the QCD improved parton model approach, there is no 
reason why the n = 1 moment of the polarized gluon distribution should not be 
considered ||21|| . In the n — *■ 1 limit, we expect in Eq.(j30|) 

A\ -> finite, A n _ -> 0, A n NS -> finite 

E\ -> 0, B n _ -> finite, £™ 5 -> . (37) 

However, ^4" , ^4^r S , and are multiplied by the following vanishing factors 

1 _ fffl\^l J 1 _ /^)\ AW% 1 f 1 _ ^ 2 (Q 2 )V- /2A+1 1 

\9 2 {p 2 )) y \ WiP 2 )) y \ \9 2 (p 2 )) j ' 

(38) 

respectively, and thus the terms proportional to P", ^4", and £>™ in Eq.(^Dj) all vanish 
in the n = 1 limit. Note that these vanishing factors are specific to the case of the 
virtual photon target, and that such factors do not appear when the target is real 
photon. Thus we find 

/ dxgj(x, Q 2 , P 2 ) = -^C n=1 + 0{a s ) . (39) 

JO 47T 2jJ 

Now due to the relation B™ = ^jBq, we have from Eq.(|34D 

C n=1 = 12ft < e 4 > (Sg + A£ G ) | n=1 . (40) 

It is noted that the combination (Bq + A™ G ) is renormalization-scheme indepen- 
dent p2| . The results in the MS scheme || ||, [L9| are 

BT 1 = 0, A-f = -2n f (41) 



The same results have been obtained by Kodaira []53 in the framework of OPE and 
RG method. He set Bq =1 = 0, observing that there is no gauge-invariant twist-2 
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gluon operator for n = 1 and obtained A™^ 1 = —2nj from the Adler-Bell-Jackiw 
anomaly. In the end, we have for the sum rule of the virtual photon structure 
function gj, 

dxgl(x,Q 2 ,P 2 ) = -— Ve 4 + 0(a s ), for A 2 < P 2 < Q 2 . (42) 
./o tt l=1 

Finally it should be pointed out that we can further pursue the QCD corrections of 
order a s to the first moment of gj [[?]] , which will be discussed in detail in [[TJJ . Here 
we just write down the result which turns out to be: 



1 dxgl{x,Q\P 2 ) = - — 

TT 



"' 1 a s (Q 2 



i=i 



TT 



2 ( g e2N2 ^ s (P 2 ) a s (Q 2 



+ 0(a 2 s ). (43) 



This result is perfectly in agreement with the one obtained by Narison, Shore and 
Veneziano in ref.@, apart from the overall sign for the definition of gj. 

The polarized virtual photon structure function gJ(x,Q 2 ,P 2 ) is recovered from 
the moments, Eq.(|30D by the inverse Mellin transformation. In Fig. 2 we have plotted, 
as an illustration, the result for n/ = 3, Q 2 = 30GeV 2 , P 2 = lGeV 2 and the QCD 
scale parameter A = 0.2GeV. The vertical axis corresponds to 

gUx,Q 2 ,P 2 )/-n f <e 4 >\n%. (44) 
tt 

Here we have shown three cases; the box (tree) diagram contribution, 

gl {B °*\x, Q 2 , P 2 ) = (2x - 1)- n, < e 4 > In % , (45) 

tt 

the LO and NLO results in QCD. We observe that the NLO QCD corrections are 
significant at large x as well as at low x. Other examples with different Q 2 and P 2 



will be reported in the forthcoming paper |15| . We have not seen any sizable change 
for the normalized structure function 

Now let us consider the real photon case P 2 = 0. The structure function can be 
decomposed as 

gl(x, Q 2 ) = gj(x, Q 2 ) | pert . + gj(x, Q 2 ) | n0 n- P ert. • (46) 
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The second term can only be computed by some non-perturbative method such as 
lattice QCD, or estimated by vector-meson-dominance model. The first term, the 
point-like piece, can be calculated in a perturbative method. Actually, it can be 
reproduced formally by setting P 2 = A 2 in Eq.(^). In Fig.3, we have plotted the 
point-like piece of gj of the real photon. The LO QCD result coincides with the 
previous result obtained by one of the authors in ||. The NLO result is qualitatively 
consistent with the analysis made by Stratmann and Vogelsang ||10|| . 

In the present paper, we have investigated polarized virtual photon structure 
function gj(x, Q 2 , P 2 ) in the NLO in QCD for the kinematical region A 2 < P 2 < Q 2 . 
The analysis has been made in the framework of the QCD improved parton model 
with the DGLAP parton evolution equations. We have shown that the behavior of 
gj(x, Q 2 , P 2 ) can be predicted up to NLO in QCD without any free parameter. In 
fact the NLO QCD corrections are significant at large x and also at small x. The 
first moment of gj for the virtual photon is non- vanishing in contrast to the case of 
real photon where the first moment vanishes. 
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Figure Caption 



Figure 1 

Deep inelastic scattering on a virtual photon in polarized e + e~ collision, 

e + e~ — > e + e~ + hadrons. The mass squared of the "probe" ("target") photon is 

-Q 2 (-P 2 ) ( A 2 < P 2 < Q 2 ). 

Figure 2 

The polarized virtual photon structure function gj(x, Q 2 , P 2 ) to the next-to- leading 
order (NLO) in units of ^ n f < e 4 > In % for Q 2 = 30GeV 2 , P 2 = lGeV 2 and the 
QCD scale parameter A = 0.2GeV with nf = 3. We also plot the leading order (LO) 
result (long-dashed line) and the box diagram contribution (short-dashed line). 

Figure 3 

The point-like piece of the polarized real photon structure function gj(x, Q 2 ) to the 
next-to-leading order (NLO) in units of ^ n f < e 4 > In ^ for Q 2 = 30GeV 2 , and 
the QCD scale parameter A = 0.2GeV with rif = 3. We also plot the leading order 
(LO) result (long-dashed line) and the box diagram contribution (short-dashed line). 
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